A dynamical approach to glassy materials 
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Typical properties of glassy materials are shown to be captured by a mean-field free- volume the- 
ory. Relaxation processes are supposed to be free-volume activated, and different entropy barriers 
are associated with density relaxation and shear motion. Free-volume time logarithmic relaxation, 
Kohlrausch-Williams- Watts, and power law viscosity result from the non-linear dynamics of spatially 
averaged quantities. The exponents associated with these phenomena are related to a single param- 
eter of the theory. The theory also accounts for coexistence of jamming transitions and non-linear 
rheology. 



I. INTRODUCTION 

The emergence of slow modes of relaxations is a ma- 
jor signature of the glassy behavior. [Q The first con- 
sequence of large equilibration time is that the proper- 
ties of a glass evolve on experimental timescales, the sys- 
tem ages. Despite the huge amount of work that has 
been devoted to glassy systems, a general understand- 
ing of slow relaxation processes and of their relation to 
macroscopic rheology is still lacking. For example, ex- 
perimental and numerical data converge to a stretched 
exponential shape of relaxation phenomena, the so-called 
Kohlrausch-Williams- Watts (KWW) relaxation, §, |, f§ 
while mode coupling theory predicts power law behav- 
ior. [0 When it comes to extend these results to non- 
linear rheology, the situation is even more difficult from 
the theoretical point of view, since the system is driven 
far from equilibrium. At an applied strain rate, power 
law viscosities or the emergence of a yield stress, are 
common features of complex fluids. || MD simulations of 
Lennard-Jones systems jfj], ^, [?J and recent experiments 
on colloidal glasses ||, ||| are also consistent with such 
properties. 

Recent theoretical approaches to non-linear rheol- 
ogy [F§ [14) have called on ideas originat- 
ing in mode coupling theory, or trap models of spin 
glasses. [Ill [l6| These approaches, however, fail to pro- 
vide a description of non-linear rheology that is consistent 
with experimental observations; || moreover, these 
models, based only on a phase-space picture of glassi- 
ness, do not lead to a clear identification of real-space 
physical mechanism that are responsible for an observed 
macroscopic behavior. Some other attempts have been 
directed towards phenomenological approaches, [jlTj but 
those are still unsatisfactory, and do not offer a micro- 
scopic picture, even heuristic, of the undergoing physical 
mechanisms. 

The current work lies at the convergence point of sev- 
eral recent theoretical ideas in various fields. Firstly, it 
is based on the idea by Liu, Nagel and coworkers, that 
some unique mechanism lies behind the profound simi- 
larities displayed by structural glassy materials, be they 
colloidal glasses, granular materials, foams,. . . Eq] Draw- 



ing on this idea, the fundemental assumption underly- 
ing the current work is that macroscopic properties of 
dense materials results from major dynamical properties 
of structural rearrangement of the contact network. 

Another important element of the current work comes 
from the study of density relaxation in granular materi- 
als. [|o[ |l[ |2| Under vertical tapping, granular materi- 
als present a time-logarithmic relaxation of the occupied 
volume; this is in strong analogy with free- volume re- 
laxation, commonly observed in glasses, and of central 
importance for aging. || However, in the case of tapped 
granular materials, density relaxation cannot result from 
thermal activation since there is no thermal bath: the 
logarithmic density relaxation has been explained by in- 
troducing an equation of motion for free-volume, which 
characterizes the internal state of the granular material. 

Finally, this work uses a theory introduced recently by 
Falk and Langer to account for elasto-plastic transition 
in amorphous solids. |2^, Q The so-called shear trans- 
formation zone (STZ) theory, originates from ideas by 
Argon, Spacpcn and others to describe creep in metallic 
alloys; [^5[ |2(| |2l, it has been primarily developed 
in the framework of fracture mechanics, and provides a 
general scheme for jamming as resulting from a struc- 
turation of the material by creation of arrangements ori- 
ented along principal directions of an applied stress. STZ 
theory calls to free- volume arguments, but in fact, free- 
volume remains a parameter and enters constants of the 
theory. 

These ideas lead me to propose a general approach to 
the rheology of dense materials, leading to a simple set of 
constitutive equations. (3(J The purpose of the this article 
is to elaborate the theoretical arguments underlying these 
equations a to detail some of their consequences. 

Two types of internal state variables are used: free- 
volume, which is an isotropic property of a molecular 
packing, and populations of arrangements, which are re- 
lated to the distortion of the contact network between 
molecules (hence introduce anisotropy). Those variables 
evolves dynamically by free-volume activated processes: 
two types of rearrangements are associated to compaction 
and shear motion. Free-volume enters activation factors 
that control both types of elementary rearrangements, 
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and it is also a dynamical quantity. 

KWW relaxation has been shown to result from the ex- 
istence of a distribution of timescales; ||, U and it seems 
that this is often thought to be a necessary requirement. 
The first contribution of the current work is to show that 
it is not necessary, and that KWW relaxation also result 
from a more simple mechanism: the interplay between 
free-volume logarithmic relaxation and free-volume acti- 
vated shear deformation. 

The second contribution of this work is to establish 
a relation between KWW relaxation, time-logarithmic 
density relaxation, and power law viscosity of a strongly 
sheared glassy material. A single parameter of the theory 
determines these three types of phenomena. Power law 
viscosity will be shown to result from the existence of 
shear induced dilatancy. I propose one of the simplest set 
of equations that can account for those three properties 
is a single framework. 

The third contribution of this work is to present a 
theoretical framework which accounts for two coexisting 
mechanisms for jamming. In mode-coupling theory, jam- 
ming results from an entropy crisis: at low temperature, 
the system supports any applied stress because no motion 
of molecules is allowed. Resistance to shear is therefore 
an intrinsic property of a material, determined by the 
thermodynamic parameters of the system, and results 
from a phase transition. In STZ theory, jamming results 
from a structuration of the contact network: jamming 
requires some amount of plastic deformation to induce 
the creation of anisotropic structures that can support 
a static stress. Those two mechanisms are shown here 
to be incorporated in a unique theory. The difference 
between hard and soft glassy materials will be identified 
with the cases when one or the other process dominates. 
The interplay between those two mechanisms lead to a 
novel interpretation of brittleness and ductility, and pro- 
vide age-dependent dynamical yield criteria. 



II. CONSTITUTIVE EQUATIONS 
A. Preliminary settings 

The main assumption made in this work is that the 
out-of-equilibrium dynamics of a glass can be modeled in 
the framework of a mean-field approximation. By mean- 
field, I mean that spatially averaged quantities account 
for the local state of the material. Constitutive equa- 
tions are proposed that relate the average shear stress, 
to the average shear strain. The internal state of the ma- 
terial is determined by free-volume and by densities of 
arrangements (to be defined further). 

The forces in the material are supposed to be charac- 
terized by a stress tensor, which is written 



deformation tensor is defined accordingly, and composed 
of shear, and of isotropic deformations, directly related 
to dilatancy. The average free-volume per molecule is 
denoted Vf, and is related to the total volume V by, 

V = Nv f + Nv Tcp , 

where w rcp is the volume per molecule of the material in 
random close packing, and N is the number of molecules. 

Under strong shear, soft glassy materials, like clays, 
foams, or granular materials, can achieve large deforma- 
tions; the evolution of the system is more similar to the 
flow of a liquid than to the deformation of a solid. In this 
work I present a theory of deformation that is expected 
to hold for large deformations of amorphous solids or soft 
glassy materials. Since there is a long-standing matter of 
debate in the mechanical community about linear versus 
non-linear deformations, it seems necessary to clar- 
ify this issue. I show here that the notion of inherent 
states [^D offers a simple physical interpretation of small 
and large deformations. 

At low temperature, a supercooled liquid or a glassy 
material evolves in a complex energy landscape which can 
be partitioned in domains of influence of local minima. 
For almost all geometrical configuration of the molecules, 
a unique inherent state is defined as the local minimum 
to which the systems relaxes if suddenly quenched to zero 
temperature, p^ ] and in the absence of forcing. When 
a material is forced, a constant stress or strain pulls the 
system through the energy landscape. Each geometrical 
configuration reached lies in the domain associated with 
an instantaneous inherent state. The time series of in- 
stantaneous inherent states provides a coarse-grained de- 
scription of the deformation of the material, and is taken 
here as the definition of irreversible, plastic, deformation. 

The instantaneous inherent state is a natural reference 
state for a given molecular configuration, and deviation 
from the instantaneous inherent state defines an inher- 
ent (shear) deformation e m . Starting from a known initial 
configuration, the total shear deformation, e = e m +7, re- 
sults both from inherent deformation, and from the flow 
of inherent states in the phase space. By definition, the 
deformation 7 is measured between two inherent states: 
it is not a state variable, since it is defined only from 
the knowledge of a configuration designed as "initial" . 
e m , however, is a state variable, since it is well defined 
for any instantaneous configuration of the system. More- 
over, e m is expected to be small, since it measures the 
deformation between two nearby configurations in phase 
space; elasticity results from a harmonic approximation 
for every local minimum: it permits to write the stress 
as given by a Hooke law, a — fj,e m . The equation that 
governs a can finally be written: 



(J.(e-j) 



(1) 



-P 

(7 



P 



where P is the pressure, and a the shear stress. The 



The rate of plastic deformation 7 will come out of a sta- 
tistical analysis of hopping motion in the phase-space. It 
will be written as a function of the state variables of the 
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system, <r, Vf, an others, to define the constitutive law of 
the material. 



B. Assumptions 

1. Free-volume 

This work relies on the idea that, in dense materials, 
macroscopic deformation results from free-volume acti- 
vated rearrangements. This is expected because temper- 
ature is so low that thermal activation becomes irrele- 
vant, and activated processes are controlled by entropic 
fluctuations of the free- volume. A rearrangements occurs 
if sufficient space exists in the neighborhood of a given 
point; the transition probabilities are thus determined by 
the size distribution of voids. In order to work at a mean- 
field level, the volume distribution is characterized by the 
single scalar value Vt] in order to write activation factors, 
it is necessary to provide an Ansatz for the distribution 
of voids in the material. 

In the original works by Cohen and Turnbull on free- 
volume theory, it was argued that the voids in a mate- 
rial are given by a Poisson distribution. jl3| [O], |]| In 
their picture, the motion of a molecule in a glassy mate- 
rial is permitted only if it can move into a hole of size, 
say v\. The probability to find a hole larger than vi, 
is proportional to exp(— Vi/vt): this leads to estimate 
a diffusion constant D oc exp(— vi/vf), which is a basis 
of free- volume approaches. Important differences, how- 
ever, exist between various uses and interpretations of 
free- volume theory. 

Falk and Langer use free-volume activation factors in 
the definition of STZ theory. Their argument bor- 
rowed from original ideas by Edwards and coworkers for 
granular materials: |^], Q the only extensive vari- 
able is supposed to be the volume V; the number of 
states available to the system is roughly proportional to 
(vf /h) N , where h is an arbitrary constant with dimension 
of a volume; the entropy is defined as, 



S(V, N) ~ N In 



N In 



V - JVv' cp 
Nh 



and the intensive quantity, analogous to temperature, is 
X- 

I = ^ - J_ 

X ~ dV v f 

Free-volume thus enters activation factors of the form 
exp(-«i/«/). 

In STZ theory, this discussion turns out to be some- 
what formal since free- volume is taken constant: Vf is 
absorbed in the phenomenological constants of the result- 
ing equations. It also leads to complicated expressions for 
transition rates as functions of the applied stress (because 
vi is supposed to be a function of a). 



However, free-volume should vary when the material 
dilates or contracts. This idea has been considered re- 
cently by several authors in the analysis of density re- 
laxation of granular materials subjected to vertical tap- 
ping. [^0[ In these works, free-volume is under- 
stood as a purely dynamical quantity, and relaxes by the 
motion of single grains in holes created by volume fluctu- 
ations. A Poisson distribution of volume fluctuations is 
assumed at all times: this leads to an equation of motion 
of the form, 



v f 



-Ri exp 



Vf 



(2) 



The activation factor exp [— V\/vf] is the probability that 
a hole of size larger than v± is present at a given point in 
the material; i?i is determined by the tapping frequency. 

The current work borrows from those different lines 
of thought. Drawing on the works initiated on granular 
compaction, free-volume is a dynamical quantity; voids 
are expected to redistribute fast enough in the material; 
fast means, faster than free- volume activated collective 
rearrangements themselves. This is justified because re- 
distribution of free-volume requires very tiny displace- 
ments of the molecules, and voids, considered as parti- 
cles, are expected to diffuse faster than molecules them- 
selves: around a void, there is by definition an excess of 
free- volume. 

There is an important difference, however, from the 
picture used in |2(J, ^fj where compaction results from 
the motion of a single grain in a hole. Here, borrowing on 
the argument of STZ theory, compaction results from el- 
ementary rearrangements, involving several molecules at 
a mesoscopic scale. Probabilities of transitions are esti- 
mated by entropic arguments. Due to the redistribution 
of voids, the number of states available to the system is 
supposed to be proportional to (vf/h) N at all times: the 
fast redistribution of voids allows the system to realize a 
quasi-equilibrium thermodynamic ensemble, determined 
by an intensive variable, Vf. 

The existence of an intensive variable, associated to the 
fluctuations of molecular configurations, is supported by 
the observation of an effective fluctuation-dissipation the- 
orem in sheared fluids. For reasons of clarity, and to 
be specific, I reserve the word temperature for the inten- 
sive quantity associated with the fluctuations of energy, 
or for the specific kinetic energy;therefore, I will not com- 
ply with to the recent use of "effective temperature" for 
something that is neither thermodynamic nor granular 
temperature. Nevertheless, the observation of an effec- 
tive fluctuation-dissipation theorem indicates that such 
an intensive quantity exists, without identifying what 
this quantity actually is. Here, it is identified to jjj and 
is associated with entropic fluctuations of molecular con- 
figurations at a mesoscopic scale. 

The important novelty introduced in the current work 
it that Vf enters activation factors and is also a dynami- 
cal quantity. This is in weak analogy with granular tem- 
perature, which is a temperature, and evolves dynami- 
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cally. But granular temperature couples with properties 
of the material only as far as it determines a collision 
frequency. |36| 49 while free-volume enters exponential 
prefactors in transformation rates. Another important 
element of the current theory is that the dynamics of free- 
volume results both from density relaxation by activated 
rearrangements - which exists in the absence of forcing 
- and from dilatancy induced by the macroscopic shear 
flow - which occurs when the system is driven out-of- 
cquilibrium. The non-linear coupling with the mean-flow 
will be addressed later; I will first continue the descrip- 
tion of rearrangement processes. 



2. Two types of activated processes 

The typical elementary compaction process, corre- 
sponds to the system going through a saddle point in 
phase space. The rate of such transformations is con- 
trolled by free-volume fluctuations and this saddle is 
characterized by a typical activation volume v±, which 
depends on the shape of the molecules, and on the de- 
tails of microscopic interactions. The activated hopping 
of the system through this type of saddle permits density 
relaxation, and leads to an equation of motion for vj of 
the form (||). 

Other types of elementary rearrangements occur in 
the material, in particular leading to elementary shear 
motion. Before presenting a more elaborate scheme for 
macroscopic shear deformation, let me use a very min- 
imal argument: shear strain is supposed to be propor- 
tional to shear stress: 7 = Da; where D is proportional 
to a free- volume activation factor. 

The current approach relies on the remark that the 
relative motion of molecules that permit compaction and 
shear are different; different types of collective rearrange- 
ments, should require different activation volumes. In 
phase space, shear and compaction correspond to two 
types of saddle points, and those saddle points are en- 
tropic barriers of different heights. The activation vol- 
ume for an elementary shear is denoted ^o, and a priori 
differs from v±. This leads to 



7 cx exp 



vo 

Vf 



(3) 



One major purpose of this work is to study the con- 
sequences of introducing complexity, not through a non- 
trivial distribution of timescales, |37j but through the 
existence of two types of entropy barriers of unequal 
heights. The set of equations (|l]), (|2|), and (||) will be 
shown to account for logarithmic relaxation of Vf and for 
a KWW relaxation of the stress after a strain increment. 



3. Shear deformation 



rial described is a liquid (in fact a non-Newtonian liquid). 
It is however, a very rough description of shear motion. 

The theoretical description of shear deformation that 
I present is directly inspired by STZ theory. |24| The 
first assumption on which STZ theory is based is that 
elementary shear deformation does not occur just any- 
where in the material, but in zones, where a rearrange- 
ment is permitted by the local configuration of molecules. 
A shear transformation zone (STZ) is thus defined as a 
locus within a material, where an elementary shear is 
possible. Another essential assumption is that STZ's are 
two-state systems. This can be understood as follows: An 
elementary arrangement corresponds to the opening and 
closing of some contacts between molecules; once such a 
rearrangement has occurred somewhere in the material, 
these molecules cannot shear further in the same direc- 
tion, but they can shear backwards. A symmetry is thus 
introduced by shearing motion, and (at least) two types 
of arrangements can be identified that are transformed 
into one another, as pictured here: 
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Note that an STZ is not supposed to contain only four 
molecules, but it is a mesoscopic object that accounts for 
the structure of several neighboring molecules (say of or- 
der 10). To simplify the analysis, only one pair of types 
of arrangements is considered, aligned along the princi- 
pal directions of the stress tensor. Macroscopic motion 
results from the statistics of elementary rearrangements, 
and reads: 



7 = A (R-i 



(4) 



where n± denote the number density of arrangements, 
and where R± denote the rate of transition ± — * =p. 

Those rates result from activation processes deter- 
mined by free- volume, but also by force fluctuations. A 
STZ shears if there is sufficient free-volume at its loca- 
tion, and if the local bias of the force network triggers the 
shear in the appropriate direction. In this work, force and 
volume fluctuations are supposed to be uncorrelated; the 
probabilities associated with those two types of fluctua- 
tions factorize: R± = R v (vf)R±((r). From the preceding 
discussion, R v (vf) = cxp(— vo/vf). In order to estimate 
the factors R± associated with the fluctuations of the 
force network, let me first remark that R±{o~) is a posi- 
tive, increasing function of a. Moreover, an elementary 
shear is expected to be triggered by large forces. The 
distributions of static forces in glassy systems have been 
shown to share strong similarities with the force distribu- 
tion in a granular material, [|19f| in which case, large forces 
are distributed exponentially, pi, [35[ |50| This justifies an 
exponential dependency of R± on a: 



a. Shear transformation zones. Equation (|3|) ac- 
counts for a free-volume dependent viscosity: the mate- 



R± = Rq exp 



vo 
v f 



exp 



±- 
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with i?Oi the update frequency of microscopic processes. 
Note however, that the exponential form for the activa- 
tion factors R± plays no role in the following, since these 
factors will be linearized. The parameter \x measures the 
typical stress that must be overcome to trigger a rear- 
rangement. In hard-sphere systems, it is expected to be 
proportional to the pressure P. In general, p, is essen- 
tially the energy required to break a bond between two 
molecules, times the density of bonds. Here, ft will be 
taken constant, either because P fixed, or because it de- 
pends on some interaction potential. At some point, \x 
will be set to unity, which fixes the unit of forces. 

b. Dynamics of arrangements. Following |^3| , ^4| , 
the equation of motion for the populations n± is writ- 
ten, 



n± = R, T n T — R±n± +(77 (A c — A a n±) 



(5) 



The first two terms of the rhs account for the transi- 
tions between both types of STZ. The last term accounts 
for the renewal of molecular configuration by the macro- 
scopic flow: STZ constantly appear and disappear dur- 
ing plastic deformation. The renewal of configurations is 
supposed to be proportional to the work of plastic defor- 
mation, 177, which is thus a common factor of creation 
and annhilitation terms. 

c. Shear-induced dilatancy. In this work, the dy- 
namics of rearrangements is coupled to free-volume dy- 
namics. Part of this coupling results from the fac- 
tor R v (vf) that controls shear transformation, but the 
macroscopic flow is also expected to influence the relax- 
ation of free- volume. 

Free-volume dynamics is derived in analogy with the 
dynamics of the populations n±: plastic deformations 
corresponds to the evolution of the system in the phase- 
space along trajectories that do not necessarily permit 
a minimization of free-volume. In real space, plastic de- 
formation requires the constant renewal of local arrange- 
ments, and new arrangements are local configurations of 
large entropy; by definition, their density is not opti- 
mized. Therefore, the renewal of configurations by the 
shear flow leads to dilatancy. In experiments on suspen- 
sions, for example, the system is fluidized by applying a 
strong shear which rejuvenates the glassy structure, pjj 

In order to quantify the average dilatancy induced by 
the shear flow, I assume that a fraction of the work of 
plastic deformation, a 7 in used for enthalpy production, 
Pvf. The equation for Vf is: 



-Ri exp 



Vf 



(77 



(6) 



This equation is an essential element of the current the- 
ory. The non-linear coupling with the macroscopic flow 
shows up when the system is strongly driven out-of- 
cquilibrium. The competition between self-relaxation 
and shear induced dilatancy will be shown to lead to 
a power law viscosity. 



C. Equations of motion 

From (||), the equation of motion for the total density, 
riT = n + + n_, of arrangements reads, 

fiT = o"7(2 Ac - A a n T ) ■ 

The creation and destruction terms account for an en- 
tropic mixing of molecular configurations. The previous 
equation shows that this mixing leads to a convergence 
of nx to an asymptotic value, — 2A C /A a ', therefore, 
rir = rioo defines a state of maximal entropy of molecular 
configurations. For this reason, I do not expect «t to be 
different from unless in very special circumstances, 
e.g. if some level of crystallization exists in the material. 
In particular, an amorphous solid is usually produced by 
a rapid quench of a liquid at low temperature: the sys- 
tem comes from a state of high entropy, therefore, in the 
resulting glass, ut is expected to achieve its asymptotic 
value. In this work, I assume that ut — at all times. 
The dynamics of local arrangements is thus determined 
only by the bias n_ — n_|_ between the populations of STZ. 

Equations (|4|) and (|^) are written in a more suitable 
form by introducing the variable 



n+ 



A = 



and parameters, eo = Ao A c /A a , Mo = A$A C , and Eq = 
2egRo. It is also convenient to introduce n — vi/vq, 
Ei = Ri/vq and a = A v /(Pvo), and to rescale free- 
volume as, 

v f 

x = — ■ 

vo 

The complete set of constitutive equations is finally: 



7 
A 

X 



Eq exp 

7 



( sinh ( — 



— A cosh 



(l-«,£rA) 



-Ei exp 



(7) 
(8) 
(9) 



Those equations involve two state variables, A and x 
which play very different roles. A evolves only when the 
material is sheared, therefore accounts for some sort of 
structural memory induced by deformation imposed on 
the material. On the contrary, x evolves spontaneously, 
and relaxes towards 0. Since this relaxation is slow, x also 
encodes some sort of memory, but a very different type 
of memory than A, a memory which constantly evolves 
as the system ages. 

Note that the asymptotic value for the relaxation of 
X (which here is 0) may, in fact, be temperature de- 
pendent; this question is eluded in the current work: it 
seems a very reasonable assumption for materials like 
suspensions if not for glasses below the glass transition, 
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because the logarithmic slow-down of density relaxation 
around k prevents it to reach asymptotic values far be- 
low (xoc « «). The parameters involved in the theory 
should depend on pressure and temperature, and this 
dependency results from particularities of a given ma- 
terial; for example the update frequencies R\ and Rq 
are expected to be proportional to \fT (where T is the 
temperature) in a dense hard-sphere material. ]3(| Such 
questions, however, are not the focus of this study; in this 
article, the values of the parameters are constant and can 
be understood as a given thermodynamical state of the 
system. 



III. INTERLUDE 

Equations (0-^) present non-linear couplings between 
shear motion and density relaxation. Before studying 
the complete set of equations in various tests, I start 
with two simple cases, when one or the other process 
dominates. I show here that jamming transitions are 
captured by equations and free- volume relaxation 

by equation (|9|). 



A. Free-volume relaxation 

In the absence of shear stress, a = and 7 = 0; the 
system ((t}|]) reduces to the single equation: 



X 



-Ei exp 



which accounts for density relaxation. 

At short time, this process is dominated by the initial 
value, xo- 



X 



-Ei exp 



K 

Xo 



this defines a time scale for short-time free- volume relax- 
ation: 



Xo 

T x = — ex P 



Ei 



K 

IXo 



The long time behavior is estimated from the integral 
expression: 



rx(0) 


K 


/ exp 




'x(t) 


.X. 



dx — k exp 



x(t) 



At long time, the integral is dominated by the small val- 
ues of x, i.e. by the ultimate value reached, x(t). There- 
fore, 



x(t) 



log(£7i*/«) 



(10) 



Free- volume relaxes logarithmically with time: in the ab- 
sence of any forcing, the system ages. 



B. Jamming 

In order to isolate the STZ mechanism for jamming, I 
now assume that free-volume is constant, at some value 
X; the evolution of the system is governed by equa- 
tions (0) and (||) only, and a constant stress a is applied. 

These equations admit two stationary solutions: a 
jammed state for which, e = 7 = 0, and 



A = tanh | - 



a steady deformation regime, for which, 

1 



and, 



e = Eq exp 



sin! 1 ( %: J cosh ( ^ 

M/ Moo- VA 1 



The latter solution is unstable at small stresses, in which 
case, the material jams. For large stresses, the jammed 
state loose stability, and steady deformation results. The 
two solutions exchange stability at the yield stress a y , 
which is the solution of: 



tanh ( 



1 



C. Various limits 

In the following, I will consider several limiting cases 
of the general constitutive equations (0-^)- 

In particular, the exponential dependency of 7 as a 
function of a will not be studied, although this non- 
linearity might be important in some cases. [32| Ac- 
tivation factors involving a will thus be linearized; this 
is valid under the assumption that a << /2; ft is then in- 
corporated into constants by taking it to unity. It leads 
to: 



7 
A 

X 



E exp[-l/ X ] (cr - A) 

1 (1-0QO-A) 
eo 

— Ei exp [— k/x] + a (77 



(11) 
(12) 

(13) 



I will also consider an isotropic (liquid) limit, where 
the structural anisotropy, measured by A plays no role. 
This is the case either at timescales when the dynamics 
of A is too slow, 1/eo — > (and A = 0) or because the 
typical stress 1/fJ.Q vanishes. In this case, the rheology is 
determined by a single state variable: 



7 = E exp[-l/x] cr 

X = -Si exp [-k/x] +aaj 



(14) 
(15) 
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The sets of equations and (|§]l|) will be helpful 

to emphasize the roles of A versus x- Many aspects of 
glassiness will be shown to be captured by equations ( 111 ) 
and (|l5|). Jamming, however, is not captured at this 
level of simplification, and requires the existence of the 
variable A. 



D. Experimental protocol 

A word about the preparation of the sample: at time 
t = 0, the system is quenched from a highly mixed state, 
in which, A = and x is large. The sample then ages 
during a waiting time, t w without any forcing. After this 
relaxation, some test is performed on the system. 

Two main experimental procedures consist in forcing 
the stress, or forcing the deformation. If a stress, a, is 
applied, the plastic deformation j(t) is determined by 
constitutive equations, and the total deformation is, e = 
7 + cr//x. If the deformation, e(t), is forced, the evolution 
of a is given by equation ([!]), and couples to constitutive 
equations. 

If a small perturbation is applied to the sample, the 
response of the system results from first or der terms in 

(BHD- 



the linearized constitutive equations (11 



lc 



In relaxation processes, however, the linear approxima- 
tion is not always valid at long times, |l7|, |53| because, 
the activation factors can become small, possibly smaller 
than non-linear terms. In the current theory, all non- 
linear terms involve the work of plastic deformation 177. 
Therefore, the linear analysis is valid so long as this work 
is small compared to relaxation processes. 



IV. STRESS RELAXATION AFTER A STRAIN 
INCREMENT 

Stress relaxation is monitored after a small strain step. 
The step can be understood as a ramp in the deforma- 
tion, which occurs on a very small time interval around 
t w . The strain e determines the initial value of the stress, 
<r(£j) = /j£. During the relaxation, e = 0; the subsequent 
evolution of the stress is governed by a = —^7, and ac- 
counts for the plastic adaptation of the material to the 
imposed strain. 

The relaxation modulus, G(t) = a(t)/a(t w ) is the cor- 
relator associated with rearrangement processes: it mea- 
sures the correlation between the initial state of the ma- 
terial, and its state at any following time, t. Assuming 
that the diffusion of molecules is dominated by collective 
rearrangements, stress relaxation can then be compared 
with the decay of self-intermediate scattering functions, 
which measure the spatial decorrelation resulting from 



the diffusion of particles in the dense medium. 1 54 The 



long-time behavior of correlation functions is one of the 
major properties of glassy materials. [EJ (See also recent 
reviews by Debenedetti and Stillinger, p5| or Angell et 
al. 01). 



A. Isotropic limit 

In the iso trop ic limit, the constitutive equations 
are (|lj) and fllq). Stress relaxation is governed by the 
linearized equations: 



& = —fi Eq exp 
X = -E x exp 



To determine the relaxation modulus, G = a/a(t w ), the 
only relevant initial condition is the initial value of Xi 
when the strain is applied; this value is also determined 
by the initial free-volume xo at time t = and by the 
waiting time t w . Moreover, either E\ or ^Eq can be in- 
corporated in the time scale, therefore, only two effective 
parameters determine all the possible outcomes of those 
equations: k, and ixEpjE\. While k measures the rel- 
ative heights of entropy barriers, fiEo/E\ measures the 
relative update frequencies of shear and compaction pro- 
cesses. 



1. Long time relaxation 

At long time, the relaxation of x is given by (|l(]), and 
stress relaxation verifies: 



a = —/j, Eq exp 



r i Elt 

K 



The response modulus G(t) = a(t)/a(t w ) is, 

with an exponent (3 which is directly related to the ratio 
k between the heights of entropy barriers: 



and A = 



yE 





Ei 

K 



-1/re 



The validity of linear approximation is checked by con- 
sidering the time-dependency of the work of plastic de- 
formations 177; this term should be compared with the 
first term is free-volume relaxation. 



Ei exp [—k/x] 



t 



The stress relaxation, leads to the following evolution for 
the work of plastic deformations: 

cr<7 oc t~ 1/K exp [-2^] . 

If k > 1, or > 0, the relaxation of cry is dominated 
by the exponential factor; if n < 1, the stress goes to a 
constant, but t~ x l K decays faster than l/t. In all cases, 
the non-linear term decays faster than the linear term, 
which justifies the linear approximation at all times. 
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For f3 > 0, i.e. n > 1, the stress undergoes KWW 
relaxation. This property results directly from the exis- 
tence of two types of saddle points of unequal heights, 
when hopping motion is controlled by a logarithmically 
relaxing free- volume. Measured values of exponent j3 in 
glassy materials, ranging from 0.2 to 0.5 j§, Q] are consis- 
tent with, k = vi/vq £ [1-25, 2]. The ratio k is related to 
the shape of molecules, and to the details of their inter- 
actions. There is no a priori reason, from a phase space 
point of view, to assume a particular value for k. In this 
work I will systematically study both cases in order to 
draw a complete picture of the phenomenology described 
by the proposed constitutive equations. 

If « < 1, then [3 < 0, and the previous calculation 
shows that the stress assumes a non-vanishing asymp- 
totic value after relaxation. In this case, A < 0, and the 
stress relaxation towards its ultimate constant value can 
be further estimated as a power law: 



G(t) 



exp 
exp 



\0\ 



\A\ (: 

-\A\t~M] (l + \A\t 



\0\ 



Interpreting G as a correlation function, it also means 
that the material never completely decorrelates from its 
initial state. This non-ergodicity is an indication of 
the glass transition as understood in mode-coupling the- 
ory [0; it can be characterized by the asymptotic value 
G(t — > oo), the Edwards- Anderson parameter: 



G(t — > oo) ~ exp t 



\0\ 



which increases with the age of the sample. 

Such a dependency of the Edwards- Anderson param- 
eter of the age of the sample has been observed experi- 
mentally by Bonn and coworkers on a colloidal glass at 
very low concentrations. [j56| 



2. Short time response 

At short time, if free-volume can be assumed large or 
constant, the exponential term exp [— 1/x] does not vary 
very much, and stress undergoes an exponential relax- 
ation: 



G{t) = exp 



-[iE$ exp 



1 

Xo 



(t - t w ) 



This relaxation occurs on a time scale of order, 



1 



exp 



1 

Xo 



which should be compared with the time scale of free- 
volume relaxation, t x . If r CT >> t x , free- volume relaxes 
faster than stress, and the long time, free-volume con- 
trolled, stress relaxation can be observed. If r CT << 
t x , the stress relaxes exponentially before free- volume 



achieves sufficiently low values to slow down stress re- 
laxation. The crossover between exponential and KWW 
relaxation is determined by, t g ~ t x or, 



f 



exp 



" 1 " 


Xo 


K 




~-exp 




.Xo. 


.Xo. 



For a large initial free-volume, this reduces to, 

Ei ~ // E xo ; 

KWW relaxation can be observed for a larger value of Ex , 
or for sufficiently small initial value of the free-volume, 
x(t w ), (that is for sufficiently large waiting time t w ). 

In order to look at the short and long time response of 
the system, equations (|l|), (|l4|) and ( [i"5| ) are integrated 
numerically with a Runge-Kutta algorithm with adap- 
tive timestep. The system is prepared in a dilute state, 
X — XO) which is allowed to relax during time t w , in the 
absence of any forcing (er([0, t w ]) — 0). At time t w , a 
strain e = 10~ 3 is applied and produces an elastic re- 
sponse, cr(£j) = /-ie. The dynamics of free- volume and 
stress are then monitored. 

a. Crossover stress. In a first series of tests, t w is 
set to 0, and the relaxation spectra are displayed for dif- 
ferent values of the parameter E\. Parameters Eq, fi 
and a are set to 1; K — 2 and xo = 20. The results 
are displayed figure [j| For the first value E\ = 1, only 
exponential relaxation is observed because all the stress 
has been released before x entered the logarithmic relax- 
ation. The value E\ = 20 corresponds to the crossover 
f^EoXo- KWW relaxation is observed at later stages of 
the relaxation, with a non-vanishing values of G. For a 
given xoj the stress release during the exponential part 
of the relaxation is proportional to [i Eq , hence inversely 
proportional to E±. The plateau in the relaxation spec- 
trum can be observed only for sufficiently large values 
of Ei/(ixEq) so that G is still observable by the time x 
enters the logarithmic relaxation. 

These relaxation curves are very similar to those ob- 
served in MD simulations, [Q, |>5[ 57 or in experi- 
ments; |58| the general appearance of those curves de- 
pends sensibly on the parameters considered. From a 
practical point of view, it means that it depends on the 
time window that is available, numerically or experimen- 
tally. 

b. Aging The crossover between exponential relax- 
ation can also be studied by varying the initial free- 
volume. This amounts to varying the waiting time t w 
for a given x{t = 0). In order to observe a plateau in 
the relaxation spectrum, parameter E\ is set to 20, and 
parameters Eq, /j, and a are set to 1; k = 1. In order to 
obtain an exponential relaxation for short waiting times, 
and to observe a wide variety of behavior, the initial value 
of x is set to 100. The results are displayed figure ||. With 
this choice of parameter, the curves resemble experimen- 
tal data. |, ||, [6(], H |6| 

Such curves strongly depend on the available experi- 
mental window, and on the width of the crossover region 
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FIG. 1: Numerical integration of equations ([!]), ( |l4| ) and ( |lE| ) 
for a fixed strain e = 10~ 3 . Parameters are Eo = a = u = f, 
and k = 2. The initial stress is set to a(t w ) = ue = 10~ 3 , 
and xo = 20, tu, = 0. Top: 1/x as a function of log 10 (t) for 
£i = 1,20,40,80,100 from right to left. Bottom: relaxation 
spectra with E\ increasing from left to right. 



which is unobservable in a log-lin plot for smaller val- 
ues of Ei/ (ji Eq). To emphasize this point, another set 
of relaxation spectra is presented figure ^ for values of 
the parameters: Eq = a = ji = 1, E\ = 8, and k — 2. 
The initial condition is x(0) = 100 and waiting times 
are: t w = 0, 5, 10, 15, 50, 100. The relaxation spectra are 
presented in log-lin, log-log, and log-log(log) plot; the 
crossover between exponential and KWW relaxation is 
observable around values of the relaxation modulus of 
order 10~ 5 , which is featureless in a log-lin plot. 

c. Ergodic to non-ergodic transition. I now consider 
the dependency of relaxation spectra on the parameter k. 
This parameter is of primary importance in the current 
approach and is directly related to the existence of en- 
tropy barriers of unequal heights in the phase space. The 
parameters are taken to be, Eo = a = fx = 1, Ei = 10, 
the initial condition is %(0) = 100, and the waiting is 
set to t w = 10. This choice is primarily motivated to 
obtain easily observable features on the relaxation spec- 
trum. The results are displayed figure |J. When n goes 
to 1, the system presents a transition to non-ergodic be- 
havior. For k < 1, the asymptotic value of the stress is 
no longer vanishing. At the point k = 1, the long time 



FIG. 2: Numerical integration of equations (§), (0) and 
for a fixed strain e = 10~ 3 . Parameters are Eo = a = fi = 
1, Ei = 20, and k — 2. Initial conditions are obtained by 
a quench at time t = from a dilute state (cr(0) = and 
x(0) = 100) followed by density relaxation without forcing 
during time t w ; at time t w the stress is set to o(t w ) — ue = 
10 -3 . Top: 1/x as a function of log 10 (t — t w ) with t w = 
0, 5, 10, 15, 50, 100 from bottom to top. Bottom: relaxation 
spectra for increasing t w from left to right. 



relaxation of a verifies: 
a = —/iEq exp 



E 

a ~ —u, — — er 
p Eit 



whence an asymptotic power law relaxation: 

, s. -(iE /Ei 



G(t) 



t u 



B. General case 

KWW relaxation is also achieved with the complete set 
of equations (jlT- 13). At time t — 0, A(0) = 0, and no 
constraint is applied until time t w : c([0, t w ]) = 0. During 
the ramp, the deformation of the material is essentially 
elastic, no rearrangement occurs, hence, A(i+) = 0, while 
the imposed strain determines the initial value the stress: 

During the ensuing relaxation, a = —fi'j and, for a 
small deformation, A ~ 7/eo- It results that, at all times, 
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FIG. 3: Numerical integration of equations ([!]), ( |l4| ) and ( |lE| ) 
for a fixed strain e = 10~ 3 . Parameters are Eo = Ex = a = 
p = 1, and k = 1.5. Initial conditions are obtained by a 
quench at time t = from a dilute state (er(0) = and x(0) = 
10) followed by density relaxation without forcing during time 
t w ; at time t w the stress is set to o~(t w ) = fie — 10 -3 . Top: 
1/x as a function of log 10 (t — t w ) with = 0, 5, 10, 50, 100 
from bottom to top. Bottom: relaxation spectra for increas- 
ing t w from bottom to top. There is no inconsistency in the 
large logarithmic scale for log 10 (G): the scale is determined 
up to the reduced parameter A. Inset: log 10 (— log 10 (G)) as a 
function of log 10 (t — t w ), with t w increasing from top to bot- 
tom. The spectra with the longest waiting time are already 
in the KWW regime. 



the quantity a 



H eo A is constant: 
a + fie A = e Uq 



Here, Uq — fi, but this expression may vary, depending 
on the preparation of a sample. To treat initial values of 
a and A in the greatest generality, the initial value for 
the difference a — A is denoted e So (and here, So = /i). 
At long times, the dynamics of A and a verify, 



A 



1 

eo 
1 

eo 



Eo exp 



Eo 



E x t 

K 



(a -A) 

l/re 

(a -A) 



whence, 



A 



:S exp [A' (tg-t^)] 



FIG. 4: Numerical integration of equations ([[]), (|ll|) and ( |l^ ) 
for a fixed strain e = 10~ 3 . Parameters are Eo = a = fi = 1 
and Ex = 10; k varies. Initial conditions are obtained by 
a quench at time t = from a dilute state (cr(0) = and 
x(0) = 100) followed by density relaxation without forc- 
ing during time t w ; at time t w the stress is set to (j{t w ) = 
fie = 10 -3 . Top: 1/x as a function of log 10 (t — t w ) with 
k = 1, 1.2, 1.4, 1.6, 1.8, 2 from top to bottom. Bottom: relax- 
ation spectra for increasing k from top to bottom. 



with, 



1 



and A' 



fi - 



J_\ Eo (Eh 
e J P V K 



-l/re 



The evolution in time of the quantity a— A is governed by 
the same equation as the dynamics of a in the isotropic 
limit. Therefore, all the previous discussion presented in 
the isotropic limit pertains to the general case, up to a 
constant term, which is fixed by the conservation of the 
quantity a + fi eo A. The dynamics of A and a follow. 
In particular, the evolution of A reads, 



A(t) oc e 



Uo-Soexp [A' (tg-t?)] 



1 + M £ o 

and the relaxation modulus verifies, 



G{t) oc 



Uo + fieoSo exp [A' (i£ - t? 1 )] 



1 



The isotropic expression for G{t) is recovered by taking 
eo — > oo . Note that the relaxation of the stress could 
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lead to a negative value if Uq < 0, which depends on the 
preparation of the sample. 

For (3 > (i.e. k > 1), the stress undergoes a KWW 
relaxation towards a non-vanishing value, which depends 
on the properties of the material and on the preparation 
of the sample. For (3 < power law relaxation is observed 
towards a non- vanishing value of the stress. 

In the current work, two origins are therefore identified 
for the existence of a non- vanishing stress at long times: 
the creation of mechanical anisotropy of the contact net- 
work (jamming), or the existence of a ergodic to non- 
ergodic transition for the hopping motion in phase space 
(freezing). Recent observations by Bonn et al indicate 
that for a given material, jamming could not be identi- 
fied only to a loss of ergodicity, and that some amount 
of jamming could coexist with liquid-like behavior, j^, |^| 
Such features are allowed in the current theory and show 
the importance of the STZ mechanism for soft glassy ma- 
terials. 

The validity of the linear approximation is checked by 
considering the long times relaxation of the rate of plastic 
deformation: 



7«f 



-1/k 



exp 



and the work of plastic deformation decays like 7 since a 
goes to a constant; the non-linear term in equation (|l~3| ) 
is negligible at all times. 




FIG. 5: Bifurcation diagram for equation (|l6|). e is drawn 
as a function of x, an d the strain rate e is fixed. For e > e* , 
no solution exists, \ > 0, the system blows up; for e < e* the 
system admits one stable solution. 



Eliminating a between equations ([lj]) and (|l5|), yields: 



X = -Ei ex P 
which is positive iff, 



K 

L X 



a— exp 
Eo 



(16) 



RESPONSE TO A CONSTANT STRAIN 
RATE 

A. Steady plastic flow 



ae z±l 

> e * 



E Ei 

In the stationary flow, the relation between shear strain 
and volume reads, 



When a constant strain rate e is applied, the material 
is driven strongly out-of-equilibrium; the response of the 
system involves non-linear terms. I study here the sta- 
tionary states of a sheared material, in which, 7 = e. In 
the first sections, I make a very primitive stability anal- 
ysis of those stationary states, by assuming that 7 = e 
is fixed at all times, and that the stress is given by the 
constitutive relation (]ll]) or (fbi]); this allows to decouple 
equation (0). 

The existence of an unstability in the complete set of 
equations accounts for the emergence of stick-slip motion 
in boundary lubrica tion. ]3(| This question will be briefly 
discussed in section VC , but raises so many issues that it 
will be tackled thoroughly in another publication; for the 
sake of completeness, the calculation of the Hopf criterion 
is given in appendix |A|. 

1. Isotropic limit 

In the limit ^0 00 ( n o yield stress), the relation 
between shear rate and shear stress reads (from (|l4|)): 

1' 



g = — exp 
Eq 



X 



Eq Ei 



exp 



k + 1 



2* 



For future use, I introduce the critical rate of deforma- 
tion, 



Eq Ei 



a 



The bifurcation diagram for \ is presented figure ^. If 
e > e* , then x > at all times: free- volume diverges, the 
system blows up. If e < e*, the equation for \ admits 
one solution: 



1 



X = 



2 (log e — log e* 



The existence of a stationary value for the free-volume 
indicates that for any finite rate of deformation, aging 
stops. In fact, this is not completely exact if a vanishingly 
small shear rate was applied, such that the associated 
steady free- volume is so small, that it is not accessible at 
reasonable timescales. As a consequence, for very small 
7, the apparent viscosity saturates to a constant value. 
This is a common feature of complex liquids (see e.g. 0). 
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In the steady state the value of the shear stress is, 



2. General case 



with, 



K- 1 



1 



The system behaves like a power law fluid; the exponent 
n of the power law viscosity is related to the exponent (3 
of KWW relaxation by, 



(3 



2-/3 



those two exponents are directly related to the slope k of 
1/x as a function of hit. 

The assumptions underlying the current approach do 
not permit to evaluate the activation volumes associated 
with density relaxation and shear motion. However, the 
single parameter k determines two exponents associated 
with the linear and non-linear response of a strained ma- 
terial. 

The existence of such a relation constitutes a impor- 
tant result of the current theory, although the form given 
to this relation depends on non-generic features: how 
non-linear dilatancy is introduced in the equation of mo- 
tion for free-volume; how the update frequencies (here 
taken constant) are written. In a hard-sphere material, 
for example, collision frequencies depend on free-volume 
as 1/x, and the relation between a and e is expected to 
be modified by logarithmic corrections. 

Molecular dynamics simulations of sheared fluid are 
consistent with power law rheologies for n equal to 
0.1 j|, llll or 0-2 H 0] Recent measurement of power 
law rhcology have also been obtained on a colloidal glass, 
leading to values of n ranging from 0.1 to 0.35. [|[ A 
value of (3 of order 0.4 corresponds to n = 0.25, which 
seems in reasonable agreement with the currently avail- 
able data. 

The theory of structural rearrangement described here, 
is expected to present some degree of universality, and 
hold, in a given form, for wide classes of materials. What 
this model shows is that there is an intimate relation be- 
tween time-logarithmic density relaxation, KWW relax- 
ation, and power law rheology of a dense material, and 
that this relation may depend only on generic features 
for a given class of materials. Therefore the existence of 
a relation between (3 and n could be tested by consider- 
ing families of materials. The experimental or numeri- 
cal evaluation of such a relation appears as an important 
mean of characterization of the coupling between two ma- 
jor modes of deformation. 



The constitutive equations are now (|llj-^3|) 

steady state, 



In the 



A = 



1 



plugging expression (11) for the plastic deformation in 
equation (|l^) leads to the following relation between \ 
and a in the steady state: 



exp 



Ex 



1 

Mo 



Using ( p"3] ) again, e can then be written as a function of 



El 
a 



aE Cl 
E x 



- I rf 

a 



1 

Mo 



For large a or e, the system behaves as a power law fluid, 
with the same exponent n as in the isotropic limit. How- 
ever, i vanishes for a non-zero value of er, which deter- 
mines the yield stress: 

1 



Mo 

The following expression can also be obtained for e as 
a function of x'- 



v^o Mo Ei exp[-K/2x] 



y/a y/ E a exp [k/x] + ^iMo exp [1/x] 

which is an increasing function of x> and which has a 
maximum for x ~ > o° : 



V E Mo Ei 



/a \/E a a + -Ei Mo 



For e > e*, like in the isotropic case, the shear deforma- 
tion leads to a constant increase of free-volume, and the 
system blows up. 



B. 



Transients 



When a shear rate is suddenly imposed on a mate- 
rial, the transient dynamics often lead to a first rise 
of the stress before a dynamical yield stress is reached 
at which plastic deformations begin; once plastic defor- 
mation occur, the stress relaxes towards its strain-rate 
dependent asymptotic value. J3| These features are ob- 
served in sheared bulk solids p3fl but also in sheared thin 
films, [|S4| |55[ |66f granular materials j73| or in solid fric- 
tion. [ |52[ Another generic property of transient regimes 
is that the dynamical yield stress depends on the age of 
the material: during aging, a material strengthens. 
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FIG. 6: Numerical integration of equations ([!]), ( |l4| ) and ( |lE| ) 
for a fixed strain rate e = 0.9. Parameters are Eo = Ei — 
a — /j, = 1, and values of k: re = 0.8 (solid lines) and re = 1.2 
(dashed lines). No stress is applied during the initial den- 
sity relaxation of the material from, x(0) = 10. At various 
times, t w , the strain rate is suddenly applied and the ensuing 
dynamics of l/\ (top) and a (bottom) are displayed. The 
strong glass re < 1 presents larger of the dynamical yield, 
which result from the smaller values of % reached. 



These phenomena are captured by the proposed equa- 
tions. Strain softening and strengthening result essen- 
tially from the equation of motion for free-volume and 
are captured by the isotropic equations. 

Figure |^ present some examples of transient dynamics 
obtained by numerical integration of equations ([!]), ( p^ - 
|l5| ). Parameters are set to Eq = E\ = ^ = a = 1. The 
initial value for the free-volume (at time t = 0) is set to 
x(0) = 10, and the system relaxes without forcing during 
a time t w . At time t w , a constant strain rate is applied. 
To prevent the blow-up of the system, e must be smaller 
than e*, which equals 1 with the chosen parameters; for 
this reason, the strain rate is set to e = 0.9. Two different 
values of the parameter k have been used to explore the 
two important cases when k is larger or smaller than 1. 
The stress is displayed as a function of the total time t 
since preparation of the sample. 

A dynamical yield stress clearly emerges for a suffi- 
ciently aged material, and increases with time. The yield- 
ing of the material is accompanied by a sudden rise of 
the free- volume allowing transformation rates to display 
non- vanishing values. 



e~ e 

FIG. 7: Phase diagram for equations (JTJ> , (jjj) and @, 
and k < 1. The solid line denotes the curve /Xhopf below 
which stick-slip motion occurs. /Uhopf vanishes for e = e* 
(dashed line) where steady plastic regime disappears, leading 
the system to break-up. 



These transients result from a rather complex interplay 
between the dynamics of free- volume and shear transfor- 
mation. The complete study of this process may lead to 
many developments that would divert this work from the 
main questions that I would like to address, and I prefer 
to leave the discussion here. 



C. Stick-slip instability 



The stability analysis of equations (jll - 13 ) and ([j]) 
shows the existence of a Hopf bifurcation which is re- 
sponsible for oscillatory response to a given strain rate. 
This leads to some understanding of stick-slip motion in 
sheared lubricants |J, H. H ^ or in granular 
materials. fFH f72|? 1 In this work, 



I mention this in- 
stability for the sake of completeness, but it raises many 
questions, and will require a more thorough considera- 
tion. The calculation of the Hopf criterion of bifurcation 
is given in appendix 

In the isotropic limit, the steady plastic deformation is 
stable iff /i > /ihopf; with: 



Mhopf 



E x 1 



E (k + 1) 2 V^oia 



a e 



In 



a e 



Eq Ei 



For k > 1, ^h pf < 0: there is no Hopf bifurcation. For 
k < 1, /ihopf > on the interval [0, e*], and the steady 
state becomes unstable iff /i < /Xh op f- This criterion is 
represented figure 0. For driving strain rates above e* , 
the system looses stability, because there is no stationary 
solution for x- Below e*, the Hopf criterion separates a 
domain of stability, above /Xhopf where the system reaches 
steady plastic deformation, and and domain below this 
curve where stick-slip motion is achieved. 

The introduction of the variable A in the set of equa- 
tions (IT]), (IlljO) slightly modifies this picture (see ap- 
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FIG. 8: Phase diagram for equations (|l]), ( |l4^|l3|) and n < 1. 

Parameters used are, Eq = Ei = a = eo = /Uo = 1, and 

k = 0.8. The solid line denotes the Hopf criterion Hho P { in the 

isotropic limit. The dashed lines indicates the Hopf criterion 
,± 



Mhopf w ith the complete set of equations. 



pendix ^) . Two limit curves A t hopf emerge, which delimit 
the unstable domain in the plane (e, fi). These curves are 
displayed figure ^ along with the previous curves fihopi of 
the isotropic limit. The unstable domain is a subspace of 
what it is in the isotropic limit, and the system present 
a creeping zone, close to e = 0, when steady sliding is 
stable for very small values of the strain rate. 

Some examples of stick-slip motion are displayed fig- 
ure ^|, with the set of parameters used in the phase di- 
agram of figure |^: Eq = E\ = a = 1, and k = 0.8 in 
the isotropic limit (solid lines) and e = fxo = 1 for the 
complete set of equations (dashed lines). The fast bursts 
of plastic deformation are accompanied by sudden jumps 
in the free-volume. A does not play a major role for the 
definition of the dynamical yield stress, as seen from the 
almost identical values of the maximum stress; A, how- 
ever, strongly influences smaller values of the stress as it 
shifts upwards the unstable fixed point; this also modifies 
the periodicity of stick-slip motion. 



VI. RESPONSE TO A CONSTANT STRESS 

In a creep test, a material is tested by forcing a con- 
stant stress. Since a is constant, the elastic deformation 
is e = a I [i at all times, and the plastic deformation is 
determined by the constitutive equations. 



a 




200 



200 



FIG. 9: Solid lines: Numerical integration of equa- 
tions ([!]), ([L4|) and ( |l5| ) for a fixed strain rate e = 0.1. Param- 
eters are Eq = E\ = a = fx = 1, and k = 0.8. Initial value 
of the free-volume is \ — 10. The regime of steady plastic 
deformation is unstable and leads to stick-slip motion. Fast 
relaxations of the stress result from sudden dilatancy of the 
material. Dashed lines: Same for equations (Q), ( |l4| - [l^ ) with 
the same parameters and, eo = Mo = 1- 



A. Linear response 

1. Isotropic limit 

A small constant stress a is applied at time t = t w . 
The plastic deformation of the material is determined by 
equations ( p^[|l5| ), for a fixed a. At large times, the rate 
of plastic deformation reads 



7 ~ E — 



Ext 

K 



-1/k 



This lead to a power law variation of the compliance, 
J = 7 /cr: 



If k > 1, or (3 > 0, the compliance diverges at long times; 
if K < 1, or (3 < 0, the compliance saturates at long 
times: 



j(oo) = |A|*; 



-101 



the system jams. This jamming is purely entropic, and 
does not result from a structural anisotropy of the force 
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network. Note that the asymptotic value of the compli- 
ance decays with the increasing age of the sample. 

The validity of the linear approximation must be 
checked by considering the evolution of the work of plas- 
tic deformation: 

0-7 ott- 1 ^ ; 

the linear approximation is justified so long as this term 
is negligible, when compared to the first term in free- 
volume relaxation (p"5|), which decays as n/t. When k < 
1, the work a 7 decays faster than n/t, and therefore, the 
linear approximation is valid at all times; when k > 1, 
however, this approximation fails when: 



a En 



Ext 

K 



-1/k 



n/t 



or, 



t nX = K 



E 



1/k 



a En a 2 



The time t is a typical time after which the linear ap- 
proximation fails; this time diverges when a — > 0. It will 
appear clearer in the following that, when k > 1, free- 
volume presents a non- vanishing stationary value for any 
small applied stress. This asymptotic value is determined 
by the balance between density relaxation and shear in- 
duced dilatancy. The failure of the linear approximation, 
for k > 1, is due to the existence of a non-zero stationary 
value of x I0r anv applied stress. 



and A as, 



A 



E x t 

hi 



-1/k 



(a -A) 



From this equation, the relaxation of the variable a — A 
is easily obtained: 



a - Ait) ~ (a - A ) exp [A" (t£ - t?)} 



with 



A „ = E±(Ex 
(3e \ k 



-1/k 



When k > 1, cr— A undergoes a KWW relaxation towards 
0: the material relaxes towards a jammed state in which 
the stress is supported by structural anisotropy. When 
k < 1, a — A saturates to a non- vanishing value: at 
long times, a fraction of the stress is supported by the 
bias between arrangements, measured by A, but another 
fraction is supported by the loss of ergodicity resulting 
from the freezing of elementary shear processes. 
From equation A ~ 7/eo, the compliance reads, 



J(t) ^ e 1 



Ao 

<7 



(1-cxp [A" (*£-**)]) 



The deformation relaxes to a constant value, which de- 
pends on the preparation of the sample; the relaxation is 
either a KWW or a power law relaxation depending on 
k. Note that the power law variation of the compliance 
that was found in the isotropic limit is recovered when 
A" — ► 0; for finite A" , this power law behavior shows up 
at short times, when A is still far from saturation. 



2. General case 

In the general case, the constitutive equations are ( pd[ - 
|l3| ). The constant applied stress deforms the material 
elastically, and at time no plastic rearrangement has 
yet occurred. Note that the initial value of the state 
variable A depends on the preparation of the sample. In 
the experiments by Cloitre et al, for example, the sample 
is prepared by applying a strong shear stress, above the 
yield stress during a rather long time. J5IJ In this case, 
the initial value of A is expected to be non- vanishing, and 
this leads to their observation of strain recovery. For the 
sake of generality, I will consider that the initial value 
of A is non- vanishing: A(t+) = Ao. In this section, 
I will assume that this value of A is small enough so 
that the linear approximation is still valid. When plas- 
tic rearrangements occur in the material, they induce a 
structuration of the material, and the variable A starts 
to evolve as A ~ 7/eo • 

The plastic rate of deformation evolves as, 



-1/k 



B. Jamming 

In this section I start by considering the complete set 
of equations ( pdHl3| ), and in particular, I recall how the 
dynamics of A lead to jamming. Then I show how this 
leads to an effective equation for free-volume which is 
valid in the isotropic limit or in the general case. 

A constant stress a is applied. The dynamics of A is 
governed by the following equation: 



E 



exp 



(cr - A) (1-AtocrA) 



(17) 



Since the normalized free-volume \ enters this equation 
only through a positive factor, the time-dependent value 
of x does not change the stability analysis performed at 
any fixed free- volume. The bifurcation diagram for A is 
drawn figure |l(| 

The dynamics of the variable A leads to jamming for 
any stress smaller than the yield stress p3f 

1 
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a 

FIG. 10: Bifurcation for the dynamics of the variable A as 
defined by equation Jl7| ), The intersection between the two 
branches of solutions define the yield stress a y . For a < a y , 
the stable solution is the jammed state, while it is the plastic 
flow for a > o v . E3| 



a 




X 

FIG. 11: Bifurcation diagram for the free- volume dynamics 
defined by equation ([H]), for n > 1. For any value of the 
stress in the interval [a y , a*], there exist a stationary value of 
X- The applied stress stops aging, and the material undergoes 
steady plastic deformation. Above a*, free- volume diverges 
leading to the break-up of the material. 



in which case, A = a and e = 0. Above the yield stress, 
the steady shear flow becomes a stable solution, in which 
case, 



A 



1 



Mo o 



and 



e = Eq exp 



(18) 



If a < a y , the systems jams, and the second term 
in the equation of motion for the free- volume vanishes. 
Free-volume continues to relax logarithmically in time 
while the system does not deform. 

When a > a y , the system presents a non-zero plastic 
flow, and the evolution of free- volume is governed by the 
following asymptotic equation: 



K 


+ aEo exp 


1" 








(*--) 


x_ 


. X_ 


V Mo/ 



-Ei exp 

X X V Mo/ 

At the level of this equation, the jamming resulting from 
STZ's shows up only through the term l//io in equa- 
tion (|l9|), which shifts the square of the applied stress. 
The stability analysis pertains essentially to the isotropic 
limit. From this point on, I will consider the response 
to a constant stress as governed by equations ( |l8| ) and 
( |l9| ) (which amounts to neglect the influence of A during 
transients). 



with: 




and contracts otherwise. The resulting bifurcation dia- 
gram for x is drawn figure [ll] and [l3] for k > 1 and k < 1 
respectively. 



1. Soft glasses 

For k > 1, free- volume dynamics admit a stationary 
solution for any a € [<Tj,,<t*], with, 




For values of the applied stress in the interval [cr^cr*], 
the material undergoes steady plastic deformation: it is 
ductile. The shear deformation reads, 



Ei 
a 



aE Cl 
Ei 



- (V-J- 

v \ Mo 



and the strain rate behaves as a power of the stress for 
large stresses. When the applied stress varies in the in- 
terval, a £ [a y ,a*], the strain rate can display a whole 
interval of values 



C. Brittleness and ductility 

I now study the free-volume dynamics which result 
from equation (p~S|) . The material dilates iff a > a*, 



e e 



0,E { 



1 



cr 



Moct 



For larger values of a, above a*, free- volume diverges, 
leading to the break-up of the material. 
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a behavior of the form, 
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1 1 1 
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1,1, 



e = E [a 



fl <7 



logio(^) 




FIG. 12: Numerical integration of equations ( |l4| ) and ( ]l5| ) for 
constant stresses <r, and for k > 1. Parameters are Eo = Ei = 
a — /j, — 1, and k = 2; for these parameter, <j* = 1. Initial 
value of the free- volume is x — 100 and t w — 0. Top: 1/x as a 
function of time for a — 0.2, 0.4, 0.6, 0.8, 1 from top to bottom; 
for (7 = 1, and for all values above, free-volume diverges. 
Bottom: log 10 (J) as a function of time for a increasing from 
bottom to top. 



The resulting dynamics for the compliance and for the 
free- volume are presented figure |l^ for different values of 
the driving stress a. The first part of the dynamics is de- 
termined by the initial value of free-volume and the cor- 
responding timescale. As easily seen on the bifurcation 
diagram (pd|), for small stresses, an initially large free- 
volume decays, and slows down the plastic deformation. 
This results in a plateau for the dynamics of the com- 
pliance; a larger times, free-volume converges towards a 
constant value, and the dynamics of J resumes a steady 
increase. 

These curves are to compare with measurements of 
particle diffusion: the particle diffusion in a dense mate- 
rial is expected to be dominated by collective processes; 
the diffusive constant is then cx exp[— 1/x], and the aver- 
age motion of a particle in a glassy medium is essentially 
governed by the same equation as J. The features dis- 
played of figure [l^ are amazingly similar to the curves 
found in several recent experiments. |5l] , |6l| 

If free-volume diverges, just before break-up, the cur- 
rent equations indicate that the strain rate should display 



However, this expression must be taken with caution, 
since the strongly out-of-equilibrium behavior of a ma- 
terial during break-up, when \ is large, is beyond the 
scope of the current theory. In particular, assumption 
has been made that the distribution of voids is close to an 
equilibrium state. During break-up, the distribution of 
free-volume is expected to display strong heterogeneities 
leading to the nucleation of fractures. 



2. Hard glasses 

For k < 1, and for any a < a* , \ < 0: free- volume 
relaxes to 0. In this case, the applied stress does not 
stop aging, shear deformation is vanishingly small: the 
material creeps. At long time, the linear term dominates 
free-volume relaxation: the calculati on of li near response 
that has been presented in section VI A 1 is valid, and 
indicates that compliance saturates: the system jams. In 
general, part of this jamming is supported by the struc- 
turation of the material, which lowers the value of the 
stress that enters equation (|l9|); the remaining stress is 
supported by entropic freezing of the rearrangement dy- 
namics. 

The sign of free- volume relaxation can be inverted only 
for larger stresses, above cr*(x); it this case, the response 
depends sensitively on the value of a and on the state 
of the material at the time t w when stress is applied: 
a larger stress is required to break an older material. 
The stationary state is unstable: when a is large enough 
to invert the sign of free-volume dynamics, the material 
breaks up. In this case, there are two types of response: 
creep for low stress, and break-up for large values of the 
stress. The material is brittle. 

The resulting dynamics for free-volume and for the 
compliance are displayed figure |l4]. The saturation of 
the compliance is clearly seen for small stresses while the 
divergence occurs at a value at a = a* (x) • The shape of 
these curves changes drastically in a very small interval of 
a around a* (x) because of long transient to escape from 
the unstable fixed point. These features are commonly 
observed in creep tests. [|[ ||| 



3. Jamming and freezing 

Materials with k > 1 are liable to display ductile be- 
havior, but can break; those with k < 1 either creep or 
break-up, and could be classified as unconditionally brit- 
tle. The study of STZ equations one a one-dimensional 
elasto-plastic decohesion has shown the importance of 
dynamics in ductile to brittle transition. |74| The cur- 
rent analysis slightly contradicts those results, and indi- 
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X 

FIG. 13: Bifurcation diagram for the free- volume dynamics 
defined by equation (fl9|), for k < 1. For all a < a*, free- 
volume continues relaxation towards 0; shear deformation oc- 
curs at a vanishingly small rate; the material creeps. The sign 
of free-volume relaxation can be changed only by applying a 
stress at least larger than a*; the required stress depends on 
the age of the material. The steady solution is always unsta- 
ble: large stresses lead to break up. 



J 




FIG. 14: Numerical integration of equations ( |l4| ) and ( |li| ) 
for constant stresses a. Parameters are Eq — E\ — a — 
fi = 1, and k = 0.8. Initial value of the free-volume is \ — 
1 and t w = 0. For these values of the parameters, and of 
initial conditions, the critical value of a where the free- volume 
diverges is, ~ 1.105. Top: 1/x as a function of time for a — 
0.5, 1, 1.05, 1.075, 1.1, 1.15, 1.2 from top to bottom. Bottom: 
J as a function of time for a increasing from bottom to top. 



cates that two major types of materials can be identi- 
fied, which display different qualitative behavior: this is 
in agreement with the common idea that some materials 
arc ductile, some other brittle. 

In fact, for hard materials (n < 1), the existence of 
a strong entropic crisis precludes dynamics of rearrange- 
ments to come into play: intrinsic properties of the ma- 
terial completely dominate the dynamical process of de- 
formation. The existence of a dynamically driven transi- 
tion between ductile and brittle behavior is still expected 
for soft materials (k > 1) for which the dynamics or ar- 
rangements can fully play its role. The interplay between 
dynamics and intrinsic glassy properties of material, as 
displayed in the current approach, seems a promising per- 
spective. 

As far as fracture mechanics is concerned, this picture 
indicates that during the breaking of a hard glass, the 
motion of molecules in a material should be essentially 
localized around the surface of a crack, where there is 
some free- volume, whereas in the bulk, plastic deforma- 
tion can occur only in places where a very large stress 
in concentrated. This is apparently related to a theoret- 
ical approach of fracture mechanics, recently introduced 
by Brener and Spatschek, |75j where the dynamics of 
the fracture is controlled by diffusion of molecules along 
the surface. For soft glassy materials, on the contrary, 
the propagation of a fracture is expected to induce more 
easily observable plastic deformation in the bulk of the 
material. 



VII. CONCLUSION 

The current work relies on the introduction of dy- 
namical equations for the intensive quantity associated 
with entropic fluctuations of molecular configurations. 
In this work, this quantity is identified to free-volume, 
but other physical pictures could be given. For example, 
free- volume is related to the average number of contacts 
per molecule; in a foam, the deformation of the bubbles 
can be accounted for a effective free-volume that is de- 
termined by the packing of effective overlapping particles 
as those used in the simulations by Durian, fl7q, [77j even 
though the relative proportions of a fluid mixture is kept 
constant; in a colloidal suspension, the existence of elec- 
trostatic screening contributes to an effective size of the 
particle that depends on the packing structure. In fact, 
the assumption made in this work is rather general, and 
many situation can be expected to lead to a dynamically 
driven effective free-volume; this may explain the ubiq- 
uity of glassy behavior. 

The theoretical approach to glassy materials can be 
summarized as: how ergodicity breaks? or what er- 
godicity breaking is sufficient to account for slow modes 
of relaxation. The introduction of a whole distribution 
of timescales is a very general way to account for non- 
ergodicity. Here I assume that a much weaker ergodic- 
ity breaking is sufficient, and this is directly supported 
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by the observation that fluctuation-dissipation theorem 
holds in a weak (but, not so weak) form. This indicates 
that the underlying distributions of state variables should 
not (in fact) be too weird, and that an effective width 
accounts for its essential features. However, since the 
system is out-of-equilibrium, this width must, at least, 
become a dynamical quantity. The equation proposed 
appear as a very minimal first order closure of the dy- 
namics of an underlying distribution of voids. 

The set of curves displayed in figure ^ is very reminis- 
cent of the temperature dependency of relaxation spec- 
tra. H Obviously, the parameters of the theory should 
depend on temperature and other thermodynamical pa- 
rameters, but here I purposely do not intend to address 
this question. There are (at least) two reasons for that: 
firstly, because glassy behavior is observed in many sys- 
tems, in particular in suspensions or in granular mate- 
rials, where thermodynamical temperature is not a rele- 
vant parameter. Secondly, because it seems to me more 
interesting to extract all possible conclusions from the 
simple equations proposed without trying to incorporate 
them into a more specific setup. I do think that an 
effective dynamically driven free-volume should emerge 
from microscopic dynamics under minimal assumptions. 
Therefore, it seems to me important to try to get some 
fundamental understanding of the consequences of such 
assumptions for a fixed set of parameters, without com- 
plicating the discussion. The dependency of the param- 
eters of this theory on thermodynamic quantities is defi- 
nitely an important issue to address in the future, but it 
first requires to have some intuition about the outcome 
of those equations. This work is intended to provide a 
grid of possible behavior displayed by simple equations, 
and I hope it can be helpful for further developments. 

Finally, one major interest of the current work is, in 
my opinion, that is gives a unified picture where differ- 
ent mechanisms for jamming coexist. This coexistence of 
jamming with non-linear rheology has been observed ex- 
perimentally, but could not be captured by earlier model 
of glassy rheology. Moreover, it opens the way towards 
a description of glassy materials where thermodynamics 
(that determine the parameters of constitutive equations) 
and dynamics coexists. 
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APPENDIX A: HOPF BIFURCATION AT 
CONSTANT STRAIN RATE 

1. Isotropic limit 

In order to carry out the stability analysis of the sys- 
tem (Q), (Q, (|l^), it is convenient to introduce the vari- 
able, 



exp 



The system (|l|) , (|1J) , (|15|) reduces to two equations for 
cr and x, 



a = /i e — /i Eq exp 
—Ei exp 



X 



1 

X 

a Eq exp 



and the jacobian of this dynamical system reads: 



-fj,E (p -/.iE a(t)log[(f)} 2 
2E Q aa(/) (-« E X (j) K + E a a 2 <p) \og[<j)] 2 



In the steady state, the variable a verifies, 



a 2 = 



aE 



and 4> takes on the value, 



a e 

E a Ei 



For all K > 0, (f> s is a strictly increasing function of e. In 
terms of <p s (e), the eigenvalues A of the jacobian verify: 



A 2 + A ( M E 4 + Ei(k- 1) In [</> s ] 2 ) + 
[iE Ei( K + l) <^ +1 ln[0 s ] 2 = 







and the Hopf bifurcation is determined by a critical value 
of /x as a function of <j> s : 
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2. General case 

In the general case, the system (|l]), (]Tl)-|T3[) involves 
one additional variable (A), and the bifurcation analy- 
sis is somewhat more complicated. In the steady state, 
the relation between strain rate and free-volume can be 
rewritten as, 

VWJ^E 1 ^l 2 



Since this expression is an strictly increasing function of 
4>, it is possible to carry out the analysis in terms of 
4> — 4> s as an effective parameter. When the strain rate 
varies between and e*, cf> varies on the interval [0,1]. 
The values above 1 are meaningless because there is no 
solution to the equations of motion. 

The Hopf criterion can be written in the form, 



H = A t i 2 + Bn + C = 



with 
A = Ei 



'o 2 a e <t> 2 Uo (2E Q acl3 + E 1 ^ 4> K ) 
+ a e (E a k (/> + E 1 (1 + k) fi Q <f> K ) log(» 2 ) 
B = E E x <j) l+K (/i 2 (2E a ( j ) + E 1 fi 0( j ) K ) + aeo^o 

x (E a (2 k - 3) 4> + 2 E 1 (k - 1) fi 4> K ) log(0) 2 
+ a 2 t 2 (k - 1) {E ancj) + E x (1 + n) fi <j) K ) log(^) 4 
C = E x 2 (k-1) Mo 2k {Eaat + Eknap) log(0) 2 
x (jx + ae («- 1) logO) 2 ) 

The quantity TL has the opposite sign of the real part of 
the complex eigenvalues of the jacobian: the steady state 
is stable, iff TL > 0. 

For k > 1, and for any <f> E [0, 1], C is positive: the two 
values of /i, solutions of TL = have the same sign. This 
sign can change only at a point where C vanishes, that 
is <j> = Q or 1. The sign shared by these solutions on the 
interval 4> € [0, 1] can then be evaluated by taking the 
limit value of the coefficient B close to or 1: it is posi- 
tive. Therefore, the equation TL = has no real positive 
solution for n > 1 : there is no Hopf instability. It can be 
checked that there is no other instability, therefore, for 
K > 1, the system is stable. 

I am now considering the case when n < 1. To get a 
first idea of the phase diagram, in the parameter space 



{/z, </>}, the quantity TL can be evaluated on the line /i = 0: 
TC(fi = 0,0) = C. On this line, the positiveness of TL is 
then equivalent to, 



log(0) 2 > 



Mo 



ae Q (1 - k) 



which requires that: 



Mo 



ae (1 - k) 



<}> < ^ho P f(M = 0) = cxp 

For <f) > </>h opf the coefficient C is negative, hence, equa- 
tion TL — has a single positive solution. On the interval 
[0, </>hopf] this equation admits two solutions which are 
positive because B is negative close to <j> = 0: the so- 
lution Mhopf which vanishes at the point </>^ opf and the 
solution Mhopf which exists on the whole interval [0,1]. 
The unstability occurs for /i S [Mhopf ' Mhopf ] ■ At small 
shear rates, or small <f>, the criterion TL is dominated by: 



.4 
B 

C 

whence, 
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E E 2 a 2 e 2 Mo - l) 1+2k log($ 4 
£i 3 ae oA/ 2 (k-1) 2 3K log(0) 4 
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A E 2 ae K + l V 
— ~ Mo k-1 ! 

The solution Mhopf can then be identified to: 



close to = 0. The same expression was found in the 
isotropic limit, but here, the relation between (f> and the 
shear rate, e is different. The solution Mhopf can ^ e ex ~ 
panded as, 



Mhopf 



C _ MO 1 - K j 
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